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Abstract. 

This paper is concerned with the detailed behaviour of roll-waves under- 
going a low- frequency perturbation. We first derive the so-called Whitham 's 
averaged modulation equations and relate the well-posedness of this set of 
equations to the spectral stability problem in the small Floquet-number 
limit. We then fully validate such a system and in particular, we are able to 
construct solutions to the shallow water equations in the neighbourhood of 
modulated roll-waves profiles that exist for asymptotically large time. 
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1 Introduction 

Our goal is to perform a two-scale analysis of waves in shallow flows, the fast 
scale being locally generated by (unmodulated) periodic travelling waves de- 
scribed by two parameters (for instance wavenumber and a discharge rate) 
and the low scale obeying an averaged system for these parameters. In study- 
ing such modulated waves our motivation is three-fold. 

First, we are directly interested in a deeper understanding of such insta- 
bilities and, in this sense, this work is a piece of a wider program [U El [3l UHl 
[T6l[T5]. Here we describe the motion of a shallow flow down an inclined ramp 
by the evolution of (h, q), h{t, x) e M"*" being the fluid heigth and q(t, x) e M 
the averaged horizontal momentum at time t > and in place x e M. The 
evolution we consider is governed by the Saint- Venant equations 



taking into account viscosity {6 > 0, 6^^ being a Reynolds number), grav- 
ity (F > is a Froude number) and a turbulent (quadratic) friction along 
the bottom. The main (non obvious) physical flaw of this description lies 
probably in the form of the viscosity term dd^q, which should be replaced 
with 6dx{hdx(q/h)). But this restriction is purely motivated by writting 
convenience. The full paper would translate to the more physical case and 
actually even in this simplified case the only nonlinear system we solve is of 
quasilinear type. Roll-waves are then depicted as periodic travelling waves 
of system ([1]), going down with a velocity larger than sound speed. Our first 
purpose is thus to investigate the behaviour of the solutions to these shallow 
water equations that are low-frequency perturbations of roll-waves. 

The other two motivations consider roll- waves as representative of a wider 
class of periodic travelling waves. In the context of Lagrangian systems, 
Whitham explained how to derive an averaged system for the slow motion 
of the local parameters describing modulated periodic travelling waves |19] 
(see in particular Chapter 14). In the following we will call such averaged 
systems, Whitham's systems. In the extension of the theory to a wider 
class of systems, and in a more mathematical way, Serre brought to light a 
direct relation between low-Floquet small eigenvalues of the original system 
linearized about a given periodic travelling wave and hyperbolicity of the 
Whitham's system linearized about corresponding paremeters [TS]. Since 
Whitham's systems are first order partial differential systems, this is precisely 
a relation between some spectral stability of the wave one wants to modulate 
and the well-posedness of the corresponding averaged system. Our goal is not 
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only to extend this result to our situation (this is done in a straightforward 
way, see Lemma [T]) but also to go one step further by associating not only 
eigenvalues but also eigenvectors (see Lemma [3]). This is performed through 
a comparison of a spectral Fourier analysis of the averaged Whitham's system 
and a spectral Bloch analysis of the original Saint- Venant system. Note that 
this already provides us with at least a spectral validation of the Whitham's 
system. 

Beyond the obvious interest of Lemma [3] by itself, it is also a key-step 
for our third goal : to extend the work of Doelman, Sandstede, Scheel and 
Schneider |5], performed in the reaction-diffusion context, to system ([1]). 
Here we discuss only the inviscid part of |5] and postpone to further work 
the viscous and shock (in parameters) parts. This extension will validate the 
Whitham's system at a nonlinear level, in the sense that to any solution to 
the Whitham's system, close to a given wave, one will associate a family of 
(higher-order) approximate solutions, with a modulated profile coincinding 
with the Whitham's solution at the linear level (see Proposition [3]), that 
describes at high order a family of solutions to system ([T]) for asymptotically 
large time (see Theorem [T]). Such a kind of validation of averaged equations 
has been performed in other contexts [5[ |6] . We believe the main difference 
between these cases and ours lies in the fact that we handle an averaged 
system, our critical modes are not easily separated and we need a careful 
spectral preparation before being able to follow the strategy in [S]. Lemma [3] 
is precisely intended to fill this gap. 

Together with the usual first-order Whitham's system (see system (144|) ) 
we have discussed up to now, we also introduce an averaged second-order 
system (see system (l48i) ). both in the derivation and in the spectral parts 
of our paper. In doing so we intend both to illustrate the strength of the 
spectral study through Bloch- Fourier comparison (see Lemma H]), and to 
prepare further work extending the viscous part of [51 to our context. 

Our paper is organized as follows. In the first section, we set notations 
for the rest of this work. In the second one, we derive formally modulation 
systems and explain how to compute in the low-frequency regime modulated 
approximate solutions up to any order with respect to e where e^^ is the char- 
acteristic wavelength of perturbations. Although these formal approximate 
solutions are not directly related to the ones we justify in Proposition [3], their 
construction shed some light on the proof of Proposition [3l Then we perform 
our spectral analysis. Afterwards, assuming the needed hyperbolicity of the 
Whitham's system, we will then provide a mathematical justification of this 
system in the spirit of what was done in in the reaction-diffusion frame- 
work. Finally we explain what are the main flaws of this justification and 
what may be expected from a detailed study of the second-order modulation 
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system, postponed to further work. 



Acknowledgement: The authors warmly thank Arjen Doelman and Guido 
Schneider for kind enlightment about [5]. 

2 Set-up 

2.1 Existence of roll- waves 

We start recalling some properties of the set of periodic travelling-wave so- 
lutions to ([1]). 

We search for a periodic travelling wave in the form 

(h, q) {x,t) = (H, Q) (cot + kx) , 

with H, Q 1-periodic functions and u, k real numbers. Then u is the time fre- 
quency, k the wavenumber and (H, Q) should satisfy the ordinary differential 
system 

{coH + kQy = 0, ojQ' + k{^ + ^y = H-^ + Sk'Q". (2) 

Setting c = —u/k, the wave velocity, and integrating the first equation of ([2]) 
as cH — Q = q yields the second order differential equation 

«^c/f"-..(| + |l)' + «-(c-l)^0. (3) 

A first result, due to Dressier, yields the existence of inviscid roll- waves, 
which are necessarily discontinuous with Lax shocks as discontinuities. 

Proposition 1 (Dressier, 5 = 0) Let F > 2 and {k,q) e R* x R* fixed. 
Then there exists a unique wavespeed c = c*(q), given by 

c*(g) = gl(F5 

such that there is a 1-periodic solution H to ^s^o! o-nd this solution is 
unique (up to a translation). Moreover, the roll-wave can he alternatively 
parametrized by (h_^_,q) with h+ the non dimensional minimum fluid height 
satisfying 

/^+>]; + ^(l + VrT4F)=: 
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in such a way that 
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2 

condition: h_ + h+ = - — ; — . 

When (5 > is fixed, one can prove the existence of small amplitude 
continuous periodic travelling waves through a Hopf bifurcation argument. 
The existence result is even better as 5 ^ where the existence of large 
amplitude roll- waves, close to Dressler's roll- waves, is proved [U [Ti] . 

Proposition 2 Let F > 2 and q > 0. Then for any wavespeed c such that 
c*{q) < c < Ch = c*{q) + 0{V6), there is a unique k{c,q) such that there is a 
1-periodic solution H to o,nd this solution is unique up to translation. 
Moreover, for any fixed 5 sufficiently small, limc^c^ ks{c,q) = (the contin- 
uous roll-wave converges to a solitary wave). 

Alternatively, for any fixed 6 sufficiently small, roll-waves can be parametrized 
by {k,q); and lim^^o C5(^,^) = c*{q) (the continuous roll-wave converges to 
a Dressier roll-wave as 6 ^ 0). 

In order to get a full idea of the bifurcation scenario, the reader is referred 
to [131 [121 E] where it is described. Yet note that there the viscosity term 
is really non physical and does not even provides us with the right jump 
condition in the small viscosity limit. However, we believe the scenario is the 
right one anyway and we check it numerically in [H El |3] (with a Lagrangian 
formulation). 

We could work under the regime describe by Proposition [2] but we rather 
choose to take as an assumption that we will work in domains where the so- 
lutions to profile equation (|3|), identified when coinciding up to translation, 
are uniquely and smoothly parametrized by {k,q). Of course this assumption 
is stronger than the mere full-rank assumption in [18]. We develop conse- 
quences of this assumption at the linear level in the next subsection. 

2.2 Abstract set-up 

In this subsection, we set some abstract notations for the rest of the paper. 
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For the sake of simplicity, let us first rewrite ([T]) as 

dth + dr,q = 0, 



(4) 



dtq + da:G{h, q) = S {h, q) + Sd^q 
where we have denoted 

G'(/^,?) = 7^ + ^> S{h,q) = h-^^. (5) 

When looking for periodic travelling- wave solutions (h, q) to (jl]) with wavenum- 
ber k e R* , and frequency oj e R, one must find 1-periodic solutions {H, Q) 
to 

{uH + kQy = Q, uQ' + k{G{H,Q))' = S{H,Q) + 6k^Q" (6) 

and set (h,q){t,x) = {H,Q){ujt + kx). For writting convenience, from now 
on we denote 

g{H,Q;k) =5eQ" -k{G{H,Q)y + S(H,Q) (7) 

and for later use 

dhgiH,Q;k)[f] = -k{dHG{H,Q)f)' + dhS{H,Q)f (8) 

d,g{H,Q-k)[f] = 5k^f" -k{d,G{H,Q)f)' + d,S{H,Q)f (9) 

and 

dkg{H,Q;k) =S2kQ" -{G{H,Q))' . (10) 

Denote also wave speed c as c = —ujjk. Then system ([6]) may be reduced 
to : there exists ^ e R such that cH — q = q and i/ is a 1-periodic solution 
to 

kc^H' + g{H,cH-q;k) = 0. (11) 

We will work in a context where, once {k,q) e R* x R fixed, profile equa- 
tion (ITT]) possesses a 1-periodic solution for one and only one speed c{k,q) 
(therefore one frequency u{k,q) = —kc{k,q)) and, for c = c{k,q), equa- 
tion (ITTl) possesses a unique 1-periodic solution, up to translation. Func- 
tions c and u are smooth and we also make a smooth choice of corre- 
sponding solution H{k,q) = H{-]k,q). Accordingly let us also denote 
Q{k,q) = Q{-;k,q) = c{k,q)H{- ; k,q) —q. All these assumptions are of 
course generic and justified at least in the small- amplitude regime. 
Now differentiating equation (|TT|) in (A;, ^)- variables leads to 



C{k^,q,)[dH{k^,q^).{k,q)] = [du{k,,q^) ■ {k,q)] A-{k^,q^) 

+ kA''{k.,q,) + qA^{k.,q,) 
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:i2) 



or alternatively to 

+ kA^{k.,q^) + 



where q^) is the linear operator associated to the linearization of equa- 
tion f lTT]) around H = H{ki,,q^), namely 



^k^,q*)f = -^^{k*,q^)c{k^,q^)f' + dhG{H{k^,q^),Q{k^,q^);k^)[f] 
+ c{k„q^)dqg{H{k„q^),Q{k„q^);k,)[f] 



(14) 



and 



A^ik.^q^) = -2^^^^^H{k.,qJ (15) 

+ ^d,g{H{k.,q^),Q{k.,q^);k^)[H{k^,q^)] 

A^{k,,q^) = d,g{H{k.,q,),Q{k^,q,);k^)[l] (16) 

A'{k^,q^) = c{k^,q^fH{k.,qJ - dkg{Hik^,q^),Q{k.,q^);k^) (17) 



kl 



-d,g{H{k.,q,),Q{k.,q^y,k.)[H{k.,q^)] 



A'{K,q,) = A\K,q^) - c{K,qM-{K,q^) (18) 
= -c{k^,q^f H{k^,q^)' - dkg{H{k^,q^),Q{k^,q^)]k^) . 

Note that, associated to uniqueness of c once {k^,q^) fixed comes the fact 
that A'^{ki,,q^) does not belong to the range of C{ki,,q^) acting on Lp,,j,(]R), 
the space of 1-periodic functions square-integrable on ]0, 1[. 

Before going on with properties of £(fc*,g^), let us choose to denote, for 
any function /, = /(/c*,?^,). Likewise, in the modulational context, once 
fixed functions (/cqj 9o) variables {X, T), for any function /, we will denote 
by /o corresponding values ; for instance 

[Cof]{y;X,T) = [C{ko{X,T),q,{X,T))f{.;X,T)]{y). (19) 

Coming back to note that, due to translation invariance of equation 
ffTTj) . HI belongs to the kernel of Moreover is a Fredholm operator 
with index 0. Let us denote £^'^(/c*, g^) the formal adjoint operator of C^, and 
choose H'^'^{kt,,q^) in its kernel such that 

<H^{K,q^)-H:>= 1 (20) 

where < • ; • > is the scalar product on Lpgj,(]R), namely, for f,ge L'^^^ 



< f\9 > = fg 



f 

Jo 



Then, for any function / belonging to Lpgj.(R), / belongs to the range of 
on Lpgj.(R) if and only if 

<Hf;f>= 0. (21) 
And we may define in a unique way an inverse 

c-\k.,q^) : {h:y^{h:Y- 

In the following, will always be meant to be defined in this way. 

We will also need a projection on the range of This may be done in 
a natural way by defining 

^ trad. £ ^ 

Pik*^Q*)f = ^Hf;A^> ^"""^ n(^.,^J/ = f-Mm. (22) 
Then, for any / in Lpg^(M), equation 

L.h = wA:^ + / 

has a solution h e Lpg^(R) if and only ii w = —p*{f), and in this case a 
solution may be defined in a unique way hy h = >C~^n*/, and any solution 
may be written in a unique way 

h = C-^UJ + aHi 

with a e R. 

More generally note that, for any / in Lpgj,(M), equation 
C,h = kA': + qAl + wA": + / 
has a solution h e L"^^^ (M) if and only if 

w = duj^{k,s) - p^{f) 
and in this case a solution may be defined in a unique way by 

h = dH4k,q) + c:'uj , 
and any solution may be written in a unique way 

h = dH,{k,q) + C:^UJ + aHi 

with a G R. 
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Likewise, for any / in Lpgj.(R), equation 

C,h = kA^, + qAl + wA'^ + / 
has a solution h e L'^^^ (R) if and only if 

w = -kAc*ik,s) - p^{f) 
and in this case a solution may be defined in a unique way by 

h = dH,{k,q) + C^'UJ , 
and any solution may be written in a unique way 

h = dH,{k,q) + C-^UJ + aHi 

with q; G R. 

For later use, we will denote 

)C{k.,q^) = C^n. . (23) 

As we will also have to solve equations like h' = f for some / in Lpgr(R), 
it is worthwhile to note that this can be done in Lpgr(R) if and only if 

</>:= f / = 
Jo 

and a solution may be defined in a unique way by 

h{y) = f/ =: /(/)(?/) . (24) 
Jo 

Correspondingly, two more functions of {k,q) will play a major role in the 
modulation analysis: 

M{k,q) = ( H{y;k,q)dy, N{k,q) = c{k,q)M{k,q) -q . (25) 
Jo 

Note that, whereas we made some choice in the parametrization of solutions 
H{k,q), functions M and N do not depend on this choice. 

Besides the assumption of smooth parametrization by {k,q), we will al- 
most always work with {k,q) such that 

dqM{k,q) # (26) 
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which may be seen to hold both in small-amphtude and small-viscosity 
regimes. At last, let us say that we will also assume 

dkc{k,q) ^ and dgc{k,q) ^ . (27) 

The first part is motivated by the fact that sometimes we will switch to a 
(c, g)-parametrization and for this purpose we will need such an assump- 
tion but of course this assumption can be removed when we work with the 
usual (/c, g)-parametrization. The second part is intended to simplify some 
discussions when solving linearized equations, and is not crucial. 



3 Formal derivation of modulation systems 
3.1 First order Whitham's equations 



We follow the method proposed by Serre [IB] to derive Whitham's equations 
and first introduce rescaled variables {X,T) = {ex,et). This yields 

drh + dxQ = , 

^ ^ ^ nth \ S{h,q) 2 (28) 

OtQ. + OxCr{h,q) = — + eddxq ■ 

We then search for an expansion of (h, q) in the form 

iKq)=f]e\h\q^) f^^;X,r 

with (/i*, g*) 1-periodic in their first argument y eM.. Identifying 0{-) terms, 
one finds 

dT<pdyh° + dx^dyq^ = 0, 

drcpdyq'^ + dx<pdy{G{h^,q')) = S{h\q') + 5{dx<i)fdy. 

Let us denote k{X,T) = dx(j){X,T), u){X,T) = dT<PiX,T) and c = -| so 
that k is the local wavenumber, u the local frequency and c the local wave 
speed. The previous system is solved by 

c{X,T) = cik{X,T),qiX,T)), u{X,T) = u{k{X,T),q{X,T)) 

and {h^,q^)(y;X,T) = {H,Q){y; k{X,T),q{X,T)), where q{X,T) is a local 
discharge rate. Note that may be recovered from k and q if and only if 

dTk + dx{kc{k,q)) =0. (29) 
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We further identify 0(1) terms in (1251) . On the one hand, the mass 
conservation law yields 



dy{u{k,q)h' - kq') = -{drh^ + ^x?"). 



This equation has a solution if and only if 



> +dx <q > = . 



(30) 



Using the fact that chP — q^ = q, equation fISU]) can be equivalently written 

as 



where M and are defined in f l25|) . The system fl29ll3ip forms the system 
of Whitham's equations. Assumption f l2B]l is equivalent to this system being 
of evolution-type for {k,q). This is a first order differential system of partial 
differential equations: in what follows , we will study the hyperbolicity of 
such a system and relate the hyperbolicity to the stability of roll-waves in 
the small wavenumber regime. We will use this set of equations to construct 
approximate solutions to the full shallow water system in the neighbourhood 
of roll-waves on asymptotically large time intervals. 

3.2 Higher order approximations 

In the following, we show how to construct a higher-order approximation of 
solutions to fl25]) : for that purpose, we need also to expand the phase with 
respect to e just as in a classical WKB-type calculation. Solutions {h, q) to 
(128|) are then expanded in the form 



with {h^, g*) 1-periodic in their first argument ?/ e M. Identifying 0(e ^) still 
yields a differential system in the form 



dT{M{k,q)) + dx{N{k,q)) = 



(31) 




dT(j)^ dyh^ + dxr dyq" = 0, 

dyq' + dx^' dy {G{h\ g°)) = S{h\ g°) + 5{dx(t^ydy. 



As previously, we set 



ko = dx(p 



,0 



= dT<P 



,0 
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There must be some Qq such that c^hP{y) — q^{y) = Qq. This yields 

a;° = u{ko,qo), c° = c{ko,%) 

and the profile equation in the y-variable 

kocldyh^ + g{h'',coh^ -%;ko) = 

solved by h'^{y;X,T) = /cq, ^q). As already mentioned, the compatibil- 
ity condition dTdx4>^ = dxdr'f' yields an evolution equation for the local 
wavenumber k^: 

drko + dx{koCo) = Q . (32) 
Next, we identify 0[1) terms. First, we consider the mass conservation law: 

Note that the expansion of 4> with respect to e yields a new term that does 
not change the other compatibility condition 

BtMq + dxNo = 0. 

However the presence of this extra term will be necessary in order to com- 
pute an approximation to the next order of the solution to the Saint- Venant 
equations. Integrating this equation with respect to y yields 

coh" -q' = ^[dTl{Ho) + dxHQo)) + (uJi + Coki)^ + q,, (33) 

Ko Kq 

where q^ is a constant of integration depending on (X, T) which plays the 
role of a correction to the relative discharge rate ^q, / was defined in fl24p 
and 

ui = dT(p\ ki = dx<\>^. (34) 

There remains to determine For that purpose, we consider the mo- 

mentum equation to order C(l). This yields 

LqU^ = k^A\ + + wiA^ + i?o (35) 

with i?o a function that depends only on Uq, k^.q^ and is defined as 

= 5tQo + dx\G{Ho, Qo)] - ^ dx\kodyQo\ - 6 ko8y{dxQo) (36) 
cody[dTl{Ho) + dxHQo)] + ^d,g{Ho, Qo; ko)[dTl{Ho) + dxHQo)] • 

Kq 

12 



Then equation fl5Sl) has a solution if and only if 

ui - dwo(^i,9i) = - Po{Ro) (37) 

where p is the function defined in f l22l) . Once again, in order to recover (p^ 
from ki and Ui, one must impose 

drh - dx[dcooih,q,)] = -dx[po{Ro)] ■ (38) 

With our choice for ui, equation may be solved by 

= dHo{h,q,) + KoRo, (39) 



/C being the operator defined in fl23|) . We still need an equation to couple 
with f p8|) in order to determine (fci,^^): this is done by considering the mass 
conservation law to order 0{e). One finds 

dyicoh^ - q') = ^{{dTh' + dxq') + dy{. . .)). (40) 
fco 

This equation has a solution provided that 

dr <h^ > +dx < q^ >= 0. (41) 
Combining ^ and ^ yields 

8t idMo{h,q,)) + dx {dNo{h,q,)) 

= -Bt < ICqRo > -dx < CqICqRo > -dx I -r^po{Ro) 



k, -^^^^^V (42) 

+ dx {dr < HHo) > +dx < HQo) >) 

The system (138| H2|) on {ki,qi) is linear hyperbolic provided that the nonlin- 
ear system (12^ IHTI) is hyperbolic and this former hyperbolicity is discussed 
in the next section. Under the assumption of hyperbolicity, we can com- 
pute {ki,qi). Then, h^,q^,(j)^ are fully determined by {hP , q^ , (p^) . The next 
steps are done similarly: assume that we have determined (h\ q\ (p^)i<:n and 
let us compute (/i""*"^, g""*"^, 0""'"^). First, we determine g""*"^ as a function 
of /i""*"^, 0"+^, {h\q\(l)^)i<^n and a y-independent qn+i with the help of an 
equation similar to fHU]) . Indeed, this equation has a solution due to the 
compatibility condition similar to (14T]) with (/i", q^) replacing (h^, q^). Then, 
one inserts the expansion of g"""*"^ into the momentum equation written at 
order (9(e"). This yields an equation in the form 
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with Rn a function of {h\ q\ 0*)j<;„ and = dx(l)^^^, (^n+i = dT4>"'~^^- This 
equation has a solution provided that 

Un+l - da;o(^n+l,gn+l) = - Po{Rn) ■ 

This yields an evolution equation for kn+i- Then one can solve the equation 
on h"'~^^ by setting 

One finishes the construction at the (n + l)-step by writing the compatibility 
condition 

Bt < > +dx < >= 

which provides {kn+i,qn+i) with an equation similar to f H2]) for {ki,q^). As 
a consequence, we build (/c„+i, by solving a linear hyperbolic system. 
This gives a full description of an approximate solution of order 0(e^~^^) and 
completes the construction of approximate solution to (128|) up to any order 
with respect to e. 

In another section, we will rigorously validate this formal construction 
and use its principle to justify modulation system ( 12^ 



3.3 Two modulation systems 

In this subsection we extract from the analysis of previous subsections two 
nonlinear systems for local wavenumber and local discharge rate (k,q) that 
we think of great importance for the understanding of low-frequency pertur- 
bations of periodic travelling-wave solutions. 

First, equations (129| [3Ti) form the Whitham's averaged system for the 
system ([1]) : 

drk - dx[uj{k,q)] = , . 

dT[M{k,q)] + dx[N(k,q)] = ^^""^ 

which stays unchanged when coming back to physical variables (x, t) by set- 
ting {k,q){x,t) = (k,q)(eX,eTy. 

dtk — dx\uj{k,'qj\ = , . 

dt{M{k,q)] + dx{N{k,q)\ = . ^^^^ 

From (l26i) may be seen that this system is of evolution type. Note that, 
though this system does not contain second order terms, the presence of 
a viscosity in system ([T]) was proeminent in order to obtain it. Moreover, 
whereas working with a realistic viscosity term would not have changed the 
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form of tfie system (and the definition of M, and u), it may indeed have 
changed the actual values of M,N and u. 

Concerning the well-posedness of system f H3|) . the best one can expect 
is that the system is hyperbolic and thus possesses local-in-time smooth 
solutions. In the next section, we will connect this hyperbolicity to the 
spectral stabihty of periodic travelling-wave solutions to system ([T]) under 
low-frequency perturbations. Then, under this spectral stability assumption, 
we will validate system f H3|) by proving that any smooth solution to f H3|l 
describes at first order an e-family of modulated true solutions to ([T]) on 
asymptotically large time. 

To be able to go beyond the shock formation in finite time for solutions 
to fH5]) . we derive now a second-order modulation system. This will be done 

as 



by combining f l55| and f l2^ In order to do so, let us rewrite 

+ BT,q{ko,%) drgo + Bx,qiko,qo) dxQo 

where Bt and Bx are defined by 

BTik.,q,){k,q) = [2cAH, + dc. H,] {k,q) 

+ ld,g{H.,Q.;k,)[I{dH.{k,q))] - q 



(45) 



and 

Bx{k^,q*){k,q) 



= ^d,g{H,,Q.;k,)[I{H,)]dc.{k,q) 

+ [Lh, + dgG{H,,Q.)H, - 26k,Hi]dc.ik,q) 

+ [cldH, + dG{H,,Q.){dH,,cAH,) - 2ShcAH:] {k,q) 

+ f^dgg{H,,Q,;K)[I{dH^{k,q))] - 5 c,H[k 



c* + 



q 



Then, at order C(£^), the pair {k,q) 
drk - dx[^{k.,q)] = 

dT{M{k,q)] + {dxN{k,q)] = 



(46) 



+ 
+ 



= (/eo, f^o) + ^{ki; qi) satisfies 

-6dx{v{kA)[BT{kA)dT{kM) 
5dx{p{Kq){Bx{kA)dx{kM) 
-Sdx [^p{k,q)[BT{k,q)dT{k,q)] 

Sdx {^p{kMBx{k,q)dx{kM 

5dx{\<I{dH{k,q)dT{k,q)) >) 
5dx [\<I{dQ{k,q)dx{k,q)) >) 
6dT{<K{k,q)[BT{k,q)dT{k,q)] >) 

5dT{<lC{k,q)[Bx{k,q)dx{k,q)] >) 
5dx {c{k,q) < IC{k,q)[BT{k,q)dT{k,q)] >) 
6dx {c{k,q) < ICik,q)[Bxik,q)dxik,q)] >] 

(47) 
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which in physical variables turns into 



dtk - d^[uj{k,q)\ = -d^ {p{k,q)[BT{k,q)dt{k,q)]) 

- d,{p{k,q)[Bx{k,q)d,{k,q)]) 

dt[M{k,q)] + [d,N{k,q)] = -5, g)]) 

- [^p{k,q)[Bx{k,q)d.{k,q)]) 

+ d4l<I{dH{k,q)dt{k,q))>) 

+ d,{l<I{dQ{k,q)d,{k,q))>) 

- dti<}C{k,q)[BTik,q)dt{k,q)]>) 

- dt{<IC{k,q)[Bx{k,q)d,{k,q)] >) 

- d,{cik,q) < }Cik,q)[BTik,q)dt{k,q)] >) 

- 4 {c{k,q) < IC{k,q)[Bx{k,q)8:,{k,q)] >) . 

(48) 

In the next section, we will connect the issue of the well-posedness of sys- 
tem f l48p with some second-order properties of spectral stabilty of periodic 
travelling-wave solutions to system ([1]) under low-frequency perturbations. 
However we postpone nonlinear validation of this second-order modulation 
system to further work. Nevertheless in the last section we do explain what 
gain may be expected from its sudy. 

4 Spectral validation of modulation systems 

In this section, we carry out a spectral validation of systems (jH]) and f l48p 
by connecting their spectral properties to spectral properties of system ([T]) 
in the low-frequency regime. 

Let us fix {k^,, g^) e R* xM. Linearization will be performed for modulation 
systems around (A;*,g^), and correspondingly for the Saint- Venant equations 
around Moreover, as it is classical when studying the stability of 

a travelling wave, we will work in a co-moving frame, either (x — c^,t, t) or 
{k^,x + ujj,t). 

4.1 The Whitham's system and Evans function 

In this subsection, we show that the dispersion relation that determines the 
hyperbolicty of the Whitham's equations provides the principal part of the 
expansion of the Evans function, associated to the spectral stability of roll- 
waves. Namely we prove the following lemma. 

Lemma 1 There is a T ^ such that 

E(A,e") ^°=°^ TD{X,iy) + O {{\X\ + \u\f) 
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where E is the Evans function associated to the Saint- Venant equations (see 
dsn]) below) and D is the dispersion of the Whitham's system (see ([S2D below). 



For that purpose, we follow the computation carried in [L6\ and consider a 
parametrization of viscous roll- waves by (c, g) rather than {k,q). The equiv- 
alence between the two parametrizations relies on assumption (1271) . This 
assumption is fulfilled on the small- amplitude regime but degenerates in the 
small- viscosity limit. We will work right here in regimes where the assump- 
tion holds ; however note that some information may be obtained directly in 
{k, g)-variables as we will show with the Bloch analysis of the last subsection 
of this section. 

In this subsection let us then denote L{c,q) g R* a period and H( ■ ;c,q) 
a L(/c, g)-periodic function such that H{L ■ ; c,q)) is a 1-periodic solution to 
([2D with k = [L{c,q)]-\ Accordingly M{c,q) = l^^"'"'^ H {y; c,q)dy and 
more generally 



<f>i 



1 r 



f 



(49) 



Since we are working with continuous roll-waves, we may choose a para- 
metrization of H such that H(0;c,q) = M(c,q). Recall that this does not 
change the Whitham's system. This sligthly simplifies computations since 
then 

M H{L)' 



dM 



L 



dL + <<lH>L = <dH>L 



Writting (IH|) in frame (x — c+t, t) and linearizing yields 



(lL^dt{c,q) - L^d^c 
< dH^ >L^ dt(c,q) + MJ^c-d^q 



0, 
0. 



(50) 



or equivalent ly 



dL* 5i(c, q) 

< dH, >L, +—dL, 
J-ii, 



St{c,q) 



dxq 



0, 
0. 



(51) 



Searching for solutions in the form e'^*"'"'^*'^^(ci, qi), with (ci, qi) constant and 
non-zero and (A, u) e C x iR, yields dispersion relation = D(X, u) where 



D{X,iy) 



(52) 



Now we recall the construction of the Evans function given in [16]. Writ- 
ting the Saint- Venant equations in the frame (x — c^,t, t) and linearizing gives 
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a linear equation dtu — Ai,u = for m = {h,q), whose operator A^, has a 
spectrum on L°°(]R) composed of A e C such that there exist functions {h, q) 
and a e E>^ such that 



dx(q - c^h) + Xh 
[<1G{H^, Q^){h, q) - c^q] + \q 



0, 

dS{H^,Q.){h,q)+5dlq, 



(53) 



and 



(q,h,h')(L.) = a iq,h,h'){0) . (54) 

Setting Y = (q, h, h'), system fl53|) can be written as a first order differential 
system with periodic coefficients: 



Y' = A{\)Y. 



(55) 



Let \E'( • ; A) denote the fundamental resolvent matrix of (155|) . Then the Evans 
function is 

^(A,a) = det(^(L^; A) -ala). (56) 

To compute an expansion of E in the neighbourhood of (A, a) = (0, 1), we 
first choose a basis of solution to (155|) for A = and then continue this basis 
of solutions for small A. We choose Yi( • ; 0) = (c^,H[, H[, H") and Ys{-;0) = 
{Ci,dgHi, — 1, dgH^,, dqH'^). Note that Yi and Y2 are indeed independent. For 
y2( ■ ; 0); ciny choice completing the basis would do. But in order to compare 
with [TB] let us choose y2( • ; 0) = (c*/i2, /i2, ^2) where /i2( ■ ; 0) is a function 
such that /i2(-f'*-;0) belongs to the kernel of on L'^(M), provided by 
the Floquet analysis of profile equation ([3]), being associated to the Floquet 
multiplyier p < 1 given by 



p = exp^^J dG'(//*,Q* 



)(l,c.) 



These eigenvectors can be continued analytically with respect to A. 

Let us denote, for any function f, [/J^ = /(i^*) — cr/(0) and [/] = 
and perform the expansion 

fe(A),/^,(A)) = Y.\\q\M,). 

leli 

Integrating f l53|) i and performing the line substitution Li — c*L2 Li gives 



-A 



h,{X) 

Jo 



+ {a-l)q,.{X) 



l=Si=S3 
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with A°(A) = det(yi(0; A),y2(0; A),y3(0; A)) and qj{X) some constants such 
that 

C*/lj - Qj = XI{hj) + Qj . 

Note that qi{0) = ^2(0) = ^3(0) = 1- Now, expanding qj{X) = 
HieN -^'^j looking at ([53]) 2, one finds 

/:.[/i}(L,-;0)] = -C.[dM*{L.-)] - it- HM)^*[(\H.iL*-)l 

Therefore there are a, /3 such that 

hi = -dM* - it - H40))dgH. + ahl + . 

As a consequence, after column substitution Ci—X(3C2—X{q\—H^,{0))C3 
Ci, comes 



A%0)E{X,a) 



X^ J * + A(a - 1)H40) -X -X ^0 ' + a - 1 

-A[5,/f.]-((T-l)/i?(0) [hi] [d^H,] 
-A[5,//:]-(a-l)/i?'(0) M [^i/I] 

+ C(^(|A| + |a-l| 



Differentiating H{L{c, q), c, q) = H{0,c,q) and H'{L{c,q),c,q) = H'{0,c,q) 
with respect to (c, ^) leads to 

(whereas = (p - l)/i2(0;0) and = (p - l)/i2'(0; 0)). This is now a 
straightforward computation to show that for some F 7^ 

E{X,e'') = TD{X,u) + O {{\X\ + \iy\f) . 

Remark: the expansion of the Evans function found here is slightly different 
from the one derived in [16]: this latter one is not correct because of an 
uncorrect expansion of [gi(A) — c*/ii(A)] + (1 — cr)(gi — c*/ii)(0; A) with respect 
to A. 

For the sake of completeness, in an appendix, we provide a more geometric 
description of this subsection that makes it more comparable with [18] . 
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4.2 Formal asymptotics in the Whitham's system 

To make some direct formal use of Lemma [1], in this subsection we study 
formally the hyperbolicity of fH4|) in the regime of small viscosity, 5^0, that 
is when viscous roll-waves are close to Dressier roll-waves. For that purpose, 
we use g)-parametrization (as in Proposition [T]) and will suppose that 
all concerned quantities are regular with respect to {6,h_^_,q). Moreover, we 
write equations in the reference co-moving frame (x — Ci,t,t). One obtains 
asymptotically 



H,{q^) dtK + Km [jj dtq + {M^c^'{q.) - 1) d.,q = 0. 

so that (h_^.,q){x) = e''^(l,0) defines a solution for any u. Therefore is 
always a trivial eigenvalue and the 2 x 2-system is always hyperbolic. The 
Whitham's system being always hyperbolic in the vanishing viscosity limit 
6^0, one can expect, using a pertubation argument, that the Whitham's 
equations are hyperbolic for 6 > sufficiently small. This is in contrast 
with the results of Boudlal and Liapidevskii [4j that suggested that inviscid 
roll-waves are stable under long wavelength perturbation (in the sense that 
the Whitham's system is hyperbolic) only for roll-waves of limited periods. 
This discrepancy may come from the fact that the modulation procedure was 
not carried out properly in [1]. Besides, the formal asymptotics carried out 
here confirms formally the spectral analysis of inviscid roll-waves performed 
rigorously in [15j . 

Let us emphasize that trivial hyperbolicity in the formal asymptotic sys- 
tem comes from the fact that wave speed c does not depend on k (nor on 
h+) but only on q. 

In contrast, a discussion in the neighbourhood of Hopf 's bifurcation points 
using a perturbation argument from the constant case would show instabil- 
ity (see [3]). Yet, some numerical evidence shows that somewhere between 
bifurcation points and limiting homiclinics there are some stable roll-waves 
and that, in the weak stable sense of hyperbolicity of the Whitham's system, 
roll-waves are even stable up to the homoclinic travelling-waves. A detailed 
discussion of these former points (instablity close to bifurcation, numerical 
check of stability for roll- waves) may be found in [H El |3]. 

4.3 Second-order modulation and Bloch-wave analysis 



In this subsection, we validate spectral properties of f l48|) by a spectral Bloch- 
wave analysis of ([1]). Also we obtain relations between eigenvectors of the 
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Saint- Venant system and those of the modulation systems (both flUl) and 

iHD). 

We are thus lead to perform a Bloch study of the operator defined as 

A4h,q) =( rM*^'i'';'^r~i^^ 1 • (57) 

^ ^ \dhQAh]+dqg,{q\-u,dyq ) ^ ' 

Yet, the reader is referred to next section for a definition of the Bloch trans- 
form. Let us only introduce, for / e [— 7r,7r], the operator ^*(/) acting on 
Lpg^(M) and defined by 

[X(/)(/)](y) = e-'nA(e^'7(-))](^/)- 

The following lemma is a version of Lemma 2.1 in [10], stated for the 
Saint- Venant equations in Lagrangian coordinates and used to prove that, 
for viscous roll-waves, linear stability implies nonlinear stability. 

Lemma 2 The critical eigenvalues Xj,j = 1,2 of Ai, are analytic functions 
of the Floquet number I. 

The Jordan structure of the zero eigenspace of A*{0) = Ai, consists of a 
1-dimensional kernel and a single Jordan chain of height 2. The left kernel 
of A^, is spanned by the constant function (1,0) and {Hl,Q[) spans the right 
eigendirection lying at the base of the Jordan chain. 

Moreover, for \l\ sufficiently small, there exist dual right and left eigen- 
functions Wj{-,1) and Wj{-,1) of A*{1) associated with Xj{l), for j = 1,2, of 
form 

2 2 
k=l k=l 

where 

• {vk)k^i,2 and {vk)k^i,2 are dual bases of the total eigenspace of Ai,{l) 
associated with small eigenvalues, analytic in I, and such that 

7)2(-,0) = (l,0) , i;i(-,0) = (if:,Q:) ; 

• Q'~^Pj,i{l),l^j,2{l)), j = 1,2, and {ll3j^i{l), f3j^2{lj), i = 1,2 are analytic 
in I. 

The role of (1, 0) is a direct consequence of the fact that the first equation 
of the Saint- Venant system is a conservation law. As already pointed out, 
the role of {H[,Q'^) stems form translational invariance of system ([1]). The 
rest of the lemma, the analyticity issue, may be obtained in a standard way 
(see [in]) and is directly related to the existence of an averaged modulation 
system. 

We come to the main part of this subsection. 
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Lemma 3 For j = 1,2, 



( 1 

1 

V K 



gl + ildkQ* 



J 



+ 0{l) 



with {kj,qj) the first term in a low-frequency expansion of a corresponding 
eigenvector of the linearized Whitham's system. Thereby vi and V2 may be 
choosed so that 



( + ildkH, \ 

1 



Lemma 4 For j = 1, 2, 



+ 0{f), V2{;1) 



I 



A,(/) = ^,(0 + Of) 



1 +0{l). 



where iJ.j{l) is a Fourier eigenvalue of system fj48|) corresponding to fre- 
quency I. 

Lemma [3] is a consequence of the proof of Lemma HJ so we focus on the 
proof of the latter. Moreover, we choose to write this proof in a wave analysis 
spirit. Lemma m also comes with a better description of critical eigenvectors 
but we do not write it here. 

Writting fH8|) in frame {k^,x + a;*t, t) and linearizing yields, after some 
simplification in the second equation with help of the first one, 

dtk + k1dc^dj:{k,q) = - k^d^ {p^[BT,*dt(k,q)]) 

k*dx 

(^p,[k,Bx,Mk,q)]) 

kA < IidH.dtik,q)) >) + k.d, < I{k^H^) > dcMk,q)) 

- kA < kj{l) > d^q^ 

- dt (< K.{BT,MkM >) - dt (< lC.{KBx,A{Kq)\ >) 

(58) 

with 

Bx,* = Bx,* - c^Bt,* . (59) 
Looking for solutions of type e^'^'^"^ {k(i>) ,q(u)) with constant {k{i'),q(i/)) 
expanded into XiieN ^^{k'',q''), one finds at order C((|A| + \i^\Y) 



dM.dt{k,q) - M^dtk 



1 







k^ 



+ k^u 



ki,dkCi, k^dqCi, 
-1 



A:" 
if 
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From f l2Bl) stems that A = is the unique eigenvalue corresponding to z/ = 0. 
We now write 



A 



leN 



to get further information. 

With this notations, we obtain 



A° 



1 







+ 



k^,^kC^, k^,dqC^, 

-1 



k' 

T70 



Naturally we recover the same dipersion relation as for f l44j) . thus spectral 
Bloch-wave validation of fHHj) contains some spectral validation of fj44|) . 
At next order, holds 



' X^k^ + X^k^ + kAc,{k\t) 

k.p.[(^X'BT,. + Bx,.) {k',g')] 



AMM,(A;0,^0) + X''dM,{k\t) - A^^fcO - A^^A;! 

= k.{^< i{dH4k^,q^)) >) - k. (f < im > g) 

-A;, < /(I) > gO) - A;.AO (< /C.[(A05r,. + Bx,.) {k^,q')] > 



smce 



A°fc° + A;,dc,(A;°,g°) 
Moreover, using again f l60|) . one may obtain 



. 



= [-c2di/.-dG'(di/.,c.diJ0-25A:.d//:](A;0,Q0) 
-Sc^Hik^ - [^dgg{H^,Q.;k^)[I{l)] + d,G{H^,Q.)^% 

2cAH, + ^Jgg{H^,Q,; k,)[I{dH,)]\ {k'^q') 



(60) 



At 



+ AO 
- AO^o 



= Blik'^.f) + AOSi(fcO,gO) + {X'^fBlik'^.f) = 5,[A0](A;0, gO) . 

We now carry out a spectral Bloch-wave analysis of system ([1]). Writting 
the system in frame + w+t, t) and linearizing gives 



that is 



dth = k^dx{c^h - q), 

dtq = -oj^dxq + Q*; A;*)[/i] + dqg{H^,Q^; k^)[q] 

dt{h,q) = A*{h,q) . 



(61) 
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We look for solutions to f l61l) in the form e'^^~^'^^(hi^{x),q^{x)) with {hy,qy) 
1-periodic functions. We are only interested in spectrum near (A, v) = (0, 0). 
Therefore we set 

A = = ^i/'A' 

and 

First, there must be g e M such that c*/io — Qo = % and then 

C.ho = qAl = C.dgH, -qd^.A: 

which in turn imposes g = (since dqU^ ^ 0). Hence there is kP e M* such 
that 

ft* rC* 

Now there is g° e M such that 

c.h' -q^ = f + ^H. 

and then 

= k'^A'l + fAl + A° fcM^ . 

This forces 

A°fc° + kAc.{k^,t) = (62) 
and the existence of A;^ g M such that 

h' = dH4k',t) + ^h: . 

Then equation 

d,{c.h^-q^) = -{c.h^-q^) + X^h" + X^h' 

implies 

c^, < > — < q^ > = A° < /i^ > 

thus 

„0 ^Oa ^* uo\ ^0 



A" { dM,{k\7f) - j^k^ ) - if = . (63) 



Note that equations fl62f63p already provides a first-order spectral justifica- 
tion. Moreover there exists e M such that 

c.h^ -q^ = t - 9°/(l) - >^^fl{H.) 

+ X^^H, + AOgi//+ AOj(dff*(A;0,^0)) . 
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Now 
with 

This sets 

AiA;° + A°A;i + A;.dc.(A;°, g") = -Kp,[B,[X']{k',f)]. (64) 
Moreover there exists k"^ such that 

= dH,{k\t) + /C,/2^ + —H[ . 

k-k 

Now equation 

c,<h^ >- <q^ > = X° <h^ > + <h^ > 

leads to 

AMM^fcO,?") + AOdM,(A;\g^) + A° < /C,S,[A°](A;°,gO) > 

= ^'^'^^ + + t + AO < /(diJ,(A;0,gO)) > (65) 

- AO < /(i/.) > - g° </(!)> 

which completes our spectral justification. 

Indeed writting a linear system for (A;0,gO) + i/(/c^,g^), which is non-zero, 
leads to the same dispersion relation for A^ + uX^ in both cases. 

Let us comment somewhat on the spectrum we just described. For u = 
in the Bloch-wave analysis, A = is an eigenvalue corresponding to a 2 x 2- 
Jordan block, with {H[,Q[) as an eigenvector (and j-i^cH*, dcQ*) as its 
antecedent in the (c, g)-parametrization). For u e iM* small, two eigenvalues 
emerge from 0, tangent to the imaginary axis when hyperbolicity of (jH]) 
is met, with first order k^,l>X^, X^ being well-described by both modulation 
systems. Curvatures of the eigenvalue curves may then be extracted from 
system (ligj) . 

5 Nonlinear validation 

of the Whitham's system 

In this section, we prove the existence of a family of solutions to the shal- 
low water equations, close to a given roll-wave, converging towards a mod- 
ulated roll-wave profile described at first order by a solution to the inviscid 
Whitham's system. 
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Firstly we enumerate here all our assumptions. Not all of them will be re- 
called in Theorem[TJ Beyond the assumptions of smooth [k, g)-parametrization 
and that the Whitham's system is of evolution-type (see ( |26|) ) and hyperbolic, 
we will also ask for : 

• the spectrum of Ai, is of upper bounded real part ; 

• for any / e [— tt, tt] and any cut-off parameter li, the non critical part 
X(l)Vl of X(0 is invertible (see dH]), ([72]) and ([73]) for definitions). 

We write the latter points as assumptions to emphasize what we really use 
in our proof. But these former assumptions may be removed. For instance, 
the upper-boundedness is a consequence of high-frequency estimates in [IB] . 
Alternatively, invertibility assumption, as hyperbolicity, can be rigourously 
reduced to numerical investigation in an explicit finite box of eigenvalue phase 
space and then numerically checked with techniques in [H [2] [3]. 

5.1 Spaces 

Following the strategy introduced in [llj for the Ginzburg-Landau equations 
and developped in [5J for reaction-diffusion systems, we will prove the con- 
vergence to a roll-wave profile in a set of analytic functions. Indeed, in 
the hyperbolic scaling considered here, there is no smoothing effect arising 
from equations whereas, in the modulation context, some terms are neglected 
precisely because they contain more derivatives. We will correspondingly re- 
strict the class of admissible solutions to the Whitham's system and to the 
Saint- Venant equations. 

Let a > and m e M"*". We first introduce a space for solutions to the 
Saint- Venant system. Let us define 



where H'^^. denotes the classical Sobolev space of 1-periodic functions. Also, 
for any Schwartz-class function u, u is the Block transform of u defined by, 
for any (y, /) e M x [— vr, tt] 




u 



(y, /) = Ju{y, = 2 ^^^'"^^{^ + 2^^'), 



u = Tu being the Fourier transform of m, explicitely for / g M 
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Note that a justification for restricting attention to R x [— tt, tt] is that ex- 
tending definition to R x R would lead to : for any {y, /) g R x R, 

u{y,l + 27r) = e-''''yu{yj) . 

Moreover, the Bloch transform comes with an inverse formula 

u(x) = ( e*'^u(x,/)di 

and a Plancherel formula 

Ml^R) = \\u\\L-^{[-^,n],Ll,,)- 

Admissible solutions will be considered in the Banach space 

= {u : R^C'^\ue Lj{a,m)} 
endowed with norm || • j^o. defined as, for any Schwartz u, 



\Mx^ - \\u\\Lj{a,m) 



J— TT 



Due to Sobolev embedding theorems, is an algebra when equipped with 
usual multiplication provided m ^ 1. Namely, if m ^ 1, there is a C{m) > 
(independent of a) such that for any u,v e 

\\uv\\x^ ^ C{m)\\u\\x^\\v\\x^ . 



Correspondingly we introduce a space for solutions to the Whitham's 
system. First 



L^(R,C) |^(0|( 
[ Jm 

and 

y^ = {u : R^C|MGL^(a,m)} 
endowed with 11 ■ ||v» : 

II w^m 



f |«(/)|(1 + |/|)'"e'^l'ld/ . 
Jr 



Note that such u are analytic on strip {z e CW^z] < a}. Moreover, when 
m ^ 1, is also an algebra. 
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An important link between two kinds of settings is provided by the fact 
that if u is 1-periodic and the Fourier transform of v is supported in [— tt, tt] 
then 

J{uv){yJ) = u{y)T{v){l) 

and in particular, for such a low-frequency v, J{v){y,l) = J^{v){l). Con- 
versely, averaging in y yields, for a general w, 

{J{u){.^)) = F{u){l) 

for any / g R (with an extended Bloch transform). 

To make our theorem more readible we also introduce more common 
spaces : uniformly local Sobolev spaces. Let us first introduce an intermedi- 
ate space 



loc 



sup \\U\\x,x+l]\\H^(Xx,x+r\) < GO 
xeR 



endowed with norm 



U||£fm = Sn\)\U\\x,x + l\\\H'^{[x,x+l]) 

" xm. 



Then we define the subspace 



I . is continuous > . 

T ^ U\- — T) ) 



5.2 Main statement 



Before stating the main result of this section, which we will prove in the 
following, we still need to introduce some change of variables. Let {k^,,q^) e 
X M. and then write equations in usual frame (k^,x + u^,t, t). Then, for any 
phase Lp, we introduce a time-dependent change of variable X'^ defined by 

X^{y,t) = y - ip{y,t) . 

If dy(p is small enough, it can be inverted into satisfying 

X'P{Y'^{x,t),t) = X . 

Note that 

. '^^^f"} , . (67) 
I - dMY'f(x,t),t) 
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In particular, whenever dyip is small, dxY^f^x^ t) is close to l + dyLpiy^^x, t),t), 
corresponding local wavenumber is close to ^^,(1 + dyip{Y'^{x,t),t)), and cor- 
responding local frequency close to 

dtip{Y^{x, t),t)+ uj.{l + dy^iV'ix, t),t)) . 

Moreover all derivatives are written as functions of Y'''{x,t) and t. 

Now the following theorem provides us with a nonlinear justification of 
the Whitham's equations 

drk + Kdx(kc(k,q) - c^k) = , . 

dTM{k, q) + k^dx {N{k, q) - cM{k, q)) = 

here written in frame {k^,X + a;*r, T). 

Theorem 1 Let (ki,q^) e M.*j_ x R. Assume that in a neighbourhood of 
(k^,,q^) system ( 168|) is of evolution type and strictly hyperbolic. 
For any a > 0, m ^ 3 and M ^ 1, there exist positive £i, r^i, Ci and Ti such 
that, for any Tq e\0,Ti], for any solution {k,q) to (^E\i on [0,Tq\ satisfying 

sup \\{k,q){-,T) - {k^,q^)\\yg ^ 7]i 
Te[0,To] 

and for all e e]0,£i[, there exist ((A;^, ^^), r|) and (^o,e such that 
sup \\{ke,q,){;t) - {k,q){e;et)\\H^ ^Ci[e + r]^], 

ie[0,To/e] 

sup \\r'^{;t)\\Hir. ^CiT]^, 
ie[0,To/e] 

sup \(po,e{t)\ ^ Ci — 
ie[0,To/e] ^ 

where 

V = sup ||(A;,^)(-,T) - (A;„^J||h- 

T£[0,To] 

and a solution (h, q) to the Saint- Venant system (|T]) such that 

sup SUp|(/l,g)(x,t) - {Happ,QappW^{x,t),t)\ ^ Ci£^, 
tE[0,To/e\ xeR 

where 

{Happ,Qapp){y,t) = {H,Q){y,ke{y,t),q^{y,t)) +r'^{y,t), 

and My,t) = ¥>oAt) + ^X-^ - ^)dz. 

Our result is partially stated in || ■ ||/^m norm, whereas we will work 
with Fourier, Bloch, or mixed types, multipliers. To fill this gap we need 
a multiplier theorem and thus refer to Lemma 5 in p^J (also stated in [5] as 
Lemma 3.6). 
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5.3 Separation of critical modes 

Our proof of the above theorem starts simultaneously rewriting the shallow 
water equations in an appropriate form so as to separate critical modes from 
others and plugging into the equations a roll-wave ansatz that we will mod- 
ulate afterwards when looking for long wavelength approximate solutions. 
Recall that, in frame (k^,x + uj^,t, t), the Saint Venant system is written as 



(69) 



dth + k^,dx{q - c^,h) = 
dtq + k^d^ {G{h, q) - c^q) = S{h, q) + 5kldlq 

or equivalently 

dth + k^^dxiq - c^h) = 
dtq + uj^d^q - Q{h, q; k^) = 

Following [5], we introduce a roll-wave ansatz for (h,q){x,t), 
Since dyip should remain small, {h, q){X'^{y,t),t) is then well-approximated 

by 

{H,,Q.){y) + {dH^,dQ.){kJy^{y,t),q{y,t)){y) + {h,q){y,t). 

Using d^niSZl) in a chain rule differentiation turns ( 1^ into a set of equa- 
tions for {h, q) and {tp, q). We added two more unknowns and therefore should 
later add two more constraints. These constraints will perform a separation 
of low-Floquet critical modes. 

As an example, note that 

would turn into 

Q* + dhQ*[dH^{k^dyip, q)] + dgQ^,[dQ^,{k^dy(f, q)] + dk<3i,k^,dy^p + 5klQ[dyip + Qn 

(taken in (y,t)), with Q-ji at least quadratic in {ki,dyip,q). 
Setting u = {h,q), the shallow water system (1^ leads to 



[Bq + B;i^(u,k^,dy(f),q)]{k^,dt(p,dtq) + dtU 



B {k^dy(p,q) - A^,u = Ti{u,k^dy(p,q) 



(70) 
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with Ai, the hnear differential operator studied in previous sections, defined in 
flFr|) and that corresponds to the hnearisation of the Saint Venant equations 
about the steady roll- wave (H^, Q*)- Recall that 



Operator Bq is given by 



Kdy{c^h - q) 



and by the fact that B^{k,q) is 



—Q[ + dQ4dy^,q) 



J 



Kdy ((c^dif* - (lQ^){k, q)) 
^hQA'^H^k^q)] + d^g^dQ^k^q)] + dkQ.k + Sk^Q'Jyk - co.dy{dQAk,q)] 

and thus may also be written as 

^ --^H[ kldc^{k,q) - H^. kAc^,{dyk, dyq) + Kdyq ^ 
k-k 

-TtQ'. kldc^k.q) + B^'idyk, dyq) 



where B^ is some differential operator with 1-periodic coefficients. At last, 
(m, k,q) is a linear differential operator whose coefficients depend at least 
linearly on (u, k, q) and 7?. is a nonlinear operator acting at least quadratically, 
informally 

B^{u,k,q) = 0{\k\ + |g| + \u\), n{u,k,q) = 0{\k\^ + + \u\''). 



We will split equation (!70|) by projecting it on low Floquet-number critical 
modes of Therefore we first introduce a projection on critical modes of 
A* for small enough Floquet numbers. 

First recall that, for / e [—n, it], ^*(/) is the operator defined by 

[X(/)(/)](2/) = e-^'nA(e^'7(-))](2/) 

so that 

A/(-,0 = X(0/(-,0 ■ 

A positive li e]0, n] can be chosen small enough so that for any / such that 
|/| ^ h the spectrum of A*{1) in a small 0-centered ball is given by two 
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spectral curves (/), j = 1,2, studied in previous section, and may be defined 
the associated ^*(/)-invariant spectral projection 



Q^(0 = ^ f(A-X(/))-MA 



(71) 



r being the boundary of our neighbourhood of 0. 

We further choose a non- increasing cut-off function, x : ^ [0, 1] 
so that, for I e M, 

' 1 if |/| ^ 1 
if 1/1^2, 



X(0 



(72) 



and define (with a slight abuse of notation), for / e [— tt,?!], truncated pro- 
jections 



vui) = Q\i) X 



m) = ^-m) 

P^f(0 = l-^mf(0- 



(73) 



To these operator-valued Bloch-symbols, we associate operators (and use an 
obvious notation for them). These operators commute and 



(1 - V')Vl, = 0, (1 - V')Vl 
We now replace system flTOj) with 



0. 



^fs^iT + ^mf^f kjy^, q)\ {Kdtip, dtq) 

dtu - A^u = VyS^ (k^dy^p, q) + V^{}Z{u, k^,dy(p, q) 

- [ViBl + Vl,Bl{u, k^dyip,q)] {k^dtif, dtq) 

supplemented with constraints 

{1-V')u = 

and 



(74) 



(75) 



supp(j'(v?),-F(g)) 



X 



(76) 



(77) 



Obviously any solution to our new formulation of the problem does provide 
us with a solution to fITOl). 
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A straightforward consequence of ([71]) is that we only need to check that 
constraint (1761) is satisfied at t = 0. In order to check whether it is so also 
for assumption ( 177|) . we need to describe precisely Vf^B^ . Of course we also 
need this detailed description to compare first equation of system (1751) with 
the Whitham's system. 

But our spectral study of previous section (see Lemmas [2] and [3]) does 
give us the needed expansion of Q". Correspondingly, we can write eigenpro- 
jection V^^ as 

'^iMy.l) = X 1 ^ j (< ^^i(-,0.^(-,/) > vi{y,l)+ < V2{-J),u{-J) > V2{y,l)) 



where again < •, • > is the scalar product of Li^^^. 

Besides, one readily obtains that, for any pair (v2,g) such that fl77j) is 
satisfied, 

J[Bl{Kv,qm) = Km{vi{l) + 0{f)) +fl{l){d^MMl) + Oil)) (78) 
therefore (using also =< diV2{Q),vi{Q) > + < V2{Q) , diVi{Q) >) 

< z)i(/), j[B^,{K^,q)m >= ( 1 + ^r,2(o ) + ^m(o 9(0, 

(79) 



< V2{l),J[B'^,{KvA)m >= ( dM. + ilS)^ [Kdyip{l)M)] 



-i < diV2{0),Vi{0) > kjyip{l), 

with formally ||%,i(/)||, ||l^i(/)|| = 0{l) and = 0{P). 

Likewise, (1771) imphes 

< v,{i),j[B''{k,q)m >= -kidc.{Hi)M))+kAi)mM)i 

< V2{l),J[B''{k,q)]{l) >= -M.kAc.{dyk{l)XW)) + kjyqil)^ 

+i < 5;?)2(0), Vi(0) > kldc^{dyk{l),dyq{l)) 

+%(0[MO,t(0], 

with ||l^,(/)||=O(0and \\i^M\ = 0{l')- 



Now, from definition of is derived for any function f 



- ^ , , ^ , is invertible, as a result, assumption ( 1771) is seen to be 
dkM^ dgM^ J 

satisfied whenever it is at t = 0. From now on we do not repeat but always 



Since 
satisf 

assume ( 177|) is satisfied 
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We still need to relate the first equation of system fl7^ with the Whitham's 
system but it is now a straightforward task. Let us split Bj(u, k,q) into 

and define iTj by i^jfil) =< Vj{-, I), /(•, /) >. Setting k = k^,^ylf and applying 
TTi to the first line of fITS]) leads to 

[l + Bl, + n,V^^,BjJu,k,q)yt{k^^) 

= -kldc.{k,q) + Bl^{k,q) + 7rlP;^f7^(n, k,q) 
-7Tir^fBl^{u, k, q)dtq - B^^dtq . 

If li is small enough, then, as long as {k,q,u) is kept small enough, this can 
be turned into 



dt{Ki^) + k1dc^,{k,q) 



l + Bl^ + Ti,V'^,Bl(u,Kq) 



-1 



{kldc^{k,q)) 



+ [1 + B{, + 7i,V'„,,BiJu, k,q)j (Bl,{k,q) + 7r^V'^,,n{u, k,q) 
1 + Bl, + 7r,V^^,BjJu, k,q)y' [{7C,V^^,Bl-^{u, k,q) + Bl,) dtq 

and, denoting the right-hand side of the former equation by 
Fk{u,k,q) - B'^g(u,k,q)dtq - B^{k,q), 

with formally \\Bl^{u,k,q){l)\\ = 0{\l\ + \u\ + \k\ + \q\), \\B^\\ = 0{l) and 
||F,(«,A;,g)|| = 0([^x|2 + |A;|2 + |g|2), gives 

dt{K^) + BlJu,k,q)dtq 



+ kidc,{k,q) + Bt{k,q) = Fk{u,k,q) 



(80) 



therefore 



dtk + dy[Biq{u,k,q)dtq 

+ k1dc,{dyk,dyq) + dy(B^{k,q)^ = dy {Fk{u,k,q)) . 

Applying rather than tti leads to 

(dM, + Sji) {dtK dtq) + n^V^^^BlJu, k,q)dt {K^) 

= -k^M^dc^(dyk, dyq) + Kdyq + B^^i^.q) 

+ 7^2p^f7^(M, k, q) - TC2V^fBl^{u, k, q)dtq 

+ i < diV2{0),vi(0) > [dtk + k^dc^dyk, dyq)] . 



^1) 
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Using flHUl) . the equation may be written as 
fdM, + BUu,k,q)){dtk,dtq) 



_ - _ _ (82) 

+ k^MAc^{dyk, dyo) - k^dyq + B^(k, q) = Fq{u, k, q) 

with formally \\B^{u,k,q)(l)\\ = 0{\l\ + \u\ + |A;| + |g|), \\B^\\ = 0{P) and 
\\Fg{u,k,q)\\ = + + |g|2). 

Again, under smallness assumptions, one can invert 

^1 U§l-^{u,k,q){l) 
duM, + §l,{u, k,q){l) dgM, + §l^{u, k,q){l) 

Thus equations ( l8T]l82i) yield an evolution system for {k,q). As a result, 
setting V = {{k,q),u), system (I75l) may be written in short as 

dtV -AV =^{V), (83) 

where A is a lower triangular linear operator written in Fourier-Bloch vari- 
ables (Fourier in {k,q), Bloch in u) a^ 



A(0 



^'^(0 
X(/) 



with b'\l) = 0{l) a bounded operator (from to 1-periodic functions) and 
A'^{1) G A^2(C) such that Af^il) = 0{l), and Af a nonlinear operator such that 
A/'(V) = (9(|Vp). From the derivation of the critical system, the spectrum of 
A'^{1) is easily seen to be given by the values j = 1, 2, of the previously 

introduced spectral curves A-,-, j = 1,2. 

In order to emphasize the mode separation, we write V = {v'^,v^) with 
= {k,q) and = u = (h, q). Then system (183|) may be written 

dtv' - A'^v' = Af'{v',v') 



where J\f^{v^,v^), J\f^{v'^,v^) = Odv^]"^ + \v^\'^) and Ar'^ is of the following form 

Ar%v'',v') = f{v^,v')dyF{v^,v') + pAf%v^,v') 

with f{v'',v'), F{v'',v') = 0{\v''\ + Iv""]), 7V'=(^;^ t;^) = + \v'\'^) and p 

is such that p(l) = 0{l). This former fact is trivial for the part coming from 



■l-With a slight abuse of notations in the use of 
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TTi through equation (I5T1) . As for the contribution of 712, it follows from the 
following fact: 

n^ihTqm =<h{.,l)> +Omh,q)] = h{l) + Omh,q)] . 

This is the end of our preparation of the system and we can now build 
a family of approximate solutions according to the desired ansatz and then 
achieve the proof with the construction of a family of solutions close to our 
family of approximate solutions. For this purpose we are now in such a 
position that we can follow the strategy explained in [5] and therefore we 
will mostly sketch the end of the proof. 

5.4 Approximate solutions 

We now fix a > 0, m ^ 3, and Tq > small enough and consider a smooth 
solution {k,q) to system fl68l) in L'^ ([0, Tq]; J^g)- For further imphcit uses, 
note that for some C(a,m) stands 

II ■ \\x^ ^ C{a,m)\\ ■ \\ya 

and that, for any e > 0, || ■ \ya is turned into || ■ \y^a by the transformation 
/ /(£ ■) so that in particular || • \ya is invariant under this transformation 
whenever 5^1. 

In order to build approximate solutions in the long-wavelength regime 
we associate to long- wavelength profiles {y^^V^) and any e > e-residuals 
Res^(y^ 1/^) and Res^(y^ V^) through 

Resl{V\V'){X,T) = [dtv'' -A.v' -bv^ -Af%v^,v')]{f,^) 

where (v^, v^) is defined by 

v\y,t) = V\ey,et), v\y,t) = V\ey,et) . 

Obviously, for a given e > 0, the above {v'^, v'^) is a solution to fl84|) if and 
only if Res^(y^ V") and Res^(V"^ V) vanish. 

The next proposition provides us with the needed approximate solutions. 
Yet, to be able to prove it we need to understand the behaviour of operators, 
defined in Bloch variables, with respect to dilatation. For this purpose, let 
us denote Da the dilatation operator, Da{f ) = f{a ■). 

First recall the diagonalisation formula 

[Tf]ix) = -^j^^e^^^[f{l)f{.,l)]{x)dl 
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where 

[f{l)g]{x) = e-'^^Tie'^-gi-m^) ■ 
Prom this one deduces 



p,-iTD,](/)(x) 



.T{el)D,]{f{;l)){x)dl. 



Note that one can not infer form this formula a Bloch transform since the 
periodicity may be lost in the process. For instance, in the case where 
T{1) = P(l, y, dy) with P a symbol 1-periodic in in the formula appears 
P {el, ^,edy). However note that even in this case thanks to the original 
periodicity the y/e-dependency is almost harmless. 

As for Bloch-Fourier operators, let us look at an operator defined through 

Tf{i) = \\r{m-my)^y 

Jo 

(with an extended definition of the Bloch transform). Then formally 
J-([D,-im](/))(0 = J]f(l+'^) f[r-(d)[e2-^-]](y)d^ 



which, when (/) = P{1, y, dy), turns into 



dy 



Note that, when / is low- frequency, it only involves bounded j/e so that, 
when moreover P is 1-periodic in y, again oscillations are harmless. 
At last, let us consider an operator defined through 

with f"(-) 1-periodic. Then 



[D,-,TD,]{f){x) = ^j'^e" 
which, when T'\y) = P{l,y), turns into 

-i=j'^e-p(d.f)/(i)di. 



(/) m' 



Again, when / is low-frequency, oscillations (here in n/e and x/e) are harm- 
less. 

We are now in position to state the following proposition. 
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Proposition 3 Let m ^ 3, M ^ 1 and uq > 0. There exists r]i > 0, si > 
and a constant Ki such that, for any < a a^, if {k,q) is a solution to 
the Whitham system such that sup^gp^T^g] \ {k,q){-,T)\ya^ ^ rji, then, for any 
< £ < £1, there exists {V^, V^) such that 

sup ||y/(£-,r)-(A;,^)(£.,r) 
sup \\V^\e-,T) 

T6[0,To] 

sup \\V:{e-,T) 

Te[0,ro] 

sup \\Resl{V,\V:){e-,T) 

T£[0,To] 

sup \\Resl{V,\V:){e-,T) 

Te[0,To] 

To prove Proposition [31 one search for (V^'^, V^) in the form 

= + eV,^ + ■■■ e'^V^ 
= + eVl + ■■■ e^Vlj 

and obtain a hierarchy of equations for (Vj ,V■')Q^j^M■ 
Kt step j, V-" is obtained as a function of {V^^ ^*^)o^i<i and by solving 
an equation (not of evolution type) through the implicit function theorem 
thanks to the inversibility of -4.*P|.. For j 7^ this function is linear in Vf^. 
Using this expression for one obtains an evolution-type equation for in 
terms of (V^'^, Vf)o<^i^j. For j = 0, the equation for Vq shares its linearization 
around [k^^q^) with the equation for {k,q). For j ^ 0, the equation is linear 
hyperbolic. 

Note that the intricated form of the proposition is a consequence of the 
fact that the norm || • \\;^a^ badly scales. Note also that, although, up to 
the expression of V?, the proof of the proposition follows the lines of the 
formal derivation of the Whitham system, the two expansions may differ 
even though they share the same starting point (k,q) at the linear level. 

5.5 From residuals to reminders 

We now look for a family ((f^, t'|))o<£<£i of solutions to system flMl) in the 
form 

(v^,,vl){x,t) = (V,^,V:){ex,et) + e'' (r';,rl){x,t) 

with (rg,r|) uniformly bounded on [0,Ti/£] (where Ti is some fixed time 
< Ti < To). 



X, 



ale ^ Kl 



e+ sup \\{k,q){-,T) - {k^,q,)\\Hr. 

T£[0,To] 



X, 



a/e ^ Kl , 



a/. sup \\{k,q){-,T) - {k^,q^)\\H. 

Te[0,To] 
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Substituting this ansatz into fl51|) yields an equation we write as 



d^rl - A.rl = MUriy^ ' ^^^^ 

Note that a hnear term b has been put in the right-hand side of the second 
equation of system fl85l) so as to deal with a diagonal form in the left-hand 
part. As usual this family of systems are solved by a fix point argument. 
For this purpose, we report the following estimates 

\mrs,rl)hf}. ^ K, + (||r,1U. + llrlll;,. ) + e^K{R,,R^) 

valid for all < a < ^ whenever IIj^^IIa'^ ^ Rc and ^ Rs, where 

goes to zero when (rji, li) goes to zero. Recall that h is a cut-off parameter 
and rji is a size parameter for {k,q). Actually Ki also (badly) depends on li 
but this latter point is not prejudicial. 

It is crucial to note that through these estimates one undergoes a loss 
of derivatives. This is a consequence of the fact that our initial nonlinear 
change of variables has turned the semilinear Saint- Venant equations into 
quasilinear ones. This former point also explains why it would not be harder 
to deal with a more physical viscosity. Unfortunately, at this stage, one can 
not exploit any smoothing coming from the linear operator But this 

is a common fact that in the process of justifying an ansatz one usually loses 
something. Here instead of losing regularity we will choose to lose analyticity. 
We will establish estimates with a width of analyticity descreasing at a linear 
pace, so that we will work on a finite time, even shorter than Tq. 

For this we introduce a smoothing operator as follows. Let us fix > 0. 
We use the fact that the spectrum of A^, is of upper bounded real part to 
choose Kq sufficiently large so that, for any / g [— 7r,7r], the real part of 
the spectrum of A^,{l) — Kq\1\ is upper bounded by —Ko\l\. We further 
define operators /cq and S^it) through their Bloch symbols k'^Q) = Kq\1\ 
and5:(f)(/) = e(^-^'0l'l. 

Then we set (r^,T|) (t) = {Ss{t)rl{t),Ss{t)r'^{t)) for times t satisfying 



^ t < 



Estimates are thus established now in || • \\xg norms. The system for the 
evolution of {r^,rl^ is similar to the previous one, with the same kind of 
estimates for nonlinear terms, but with linear operator (A'^, A*) replaced 
with (A'^ — /cq, — /cq). The reason for this change of unknowns is that now 
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one can prove for (A'^ — A;q, — /cq) a maximal regularity result in A'^- valued 
Holder spaces, readily similar to Lemma 6.3 in [S]. This was the missing part 
to close in a classical way a fix-point iteration scheme. 

The last thing we should say is that vJo.e is then recovered by integrating 
over time Ti/e equation (IHOj) . 

6 Conclusion 

In this paper, we derived formally first-order and second-order averaged equa- 
tions for shallow water fiows that describes the dynamics of modulated roll- 
waves and provided two set of justification results. 

On the one hand, we carried out a spectral stability analysis of roll-waves 
using Bloch transform in the regime of small wavenumber perturbations. We 
first related the index of stability of viscous roll-waves with the hyperbolicity 
of the first-order Whitham's equations just as it was done by Serre for general 
viscous conservation laws [TS] or Johnson, Zumbrun and Bronski for gener- 
alized Korteweg-de Vries equations However, in both latter papers, only 
the Evans function framework was used. Here we proved the Bloch trans- 
form framework is a more natural tool, by extending such stability analysises 
into two directions : we relate not only eigenvalues but also eigenvectors and 
we relate the stability of steady solutions to the second-order Whitham's 
equations with the parabolicity of spectral curves at the origin. 

On the other hand, in the spirit of what has been done for reaction- 
diffusion equations by Doelman, Sandstede, Scheel and Schneider [5], we 
justified rigorously the inviscid Whitham's equations in the natural hyper- 
bolic scaling. More precisely, we proved that, given a solution to the inviscid 
Whitham's system, there exist solutions to the viscous shallow water equa- 
tions on asymptotically large time that are close to modulated roll-waves 
whose first-order expansion is described by our solution to the Whitham's 
system at the linear level. This justification is performed under weak stabil- 
ity assumptions : at the origin, tangency to the imaginary axis of spectral 
curves. From numerical investigations, it seems that it will enable us to ap- 
ply this nonlinear justification up to the limiting homoclinic travelling waves 
(whose spectrum yields unstablility but is tangent to to the imaginary axis 
[3]). This weak stability assumption has a counterpart in the required ana- 
lyticity of solutions. 

But probably the main fiaw of this justification is that our solution to the 
Whitham's system describes the first order of the roll-wave profile only at the 
linear level. This is a consequence of both the Y"^ change of variables and the 
hyperbolic scaling. Such an issue would not occur with a diffusive scaling. 
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Yet the counterpart would be that the size of the allowed perturbations of 
(A;,^) would not be anymore 0{1) but 0{e) (recall e > is the characteristic 
wavenumber of the modulation). 

Moreover the nonlinear justifications need some regularity so that, with 
the hyperbolic scaling, there is little hope to construct modulated roll-waves 
converging at +00 to roll-waves with different wavenumbers and local dis- 
charge rate since they would correspond to shocks in parameters (k,q). How- 
ever with a parabolic scaling on may expect to justify viscous shocks of the 
second-order modulation system as modulated roll-wave profile for solutions 
to the Saint- Venant system. 

At last a second-order justification would probably also enlarge the time 
validity of the approximation from [0,To/£] with To some time imposed by 
the equations to [0,T/£^] for any fixed time T. 

For all these reasons, a natural direction would be now to justify a second- 
order modulation in a parabolic scaling as it was also done in ^ for reaction- 
diffusion equations. 

A Whitham's dispersion in a geometric way 

In this appendix, we explain how to fill the small gap between analysis 
in subsection 14.11 and the one in [IB] (for different equations). We just 
need to replace our straigthforward computations with more geometric ones. 
Yet we keep the framework of subsection 14.11 : frame {x — c^,t,t), (c, g)- 
parametrization. . . 

We turn 2 x 2-determinant giving Whitham's dispersion into a 3 x 3- 
determinant. This is done interpreting the set of periodic travelling- wave 
solutions submanifold. 

The tangent space of the set of periodic travelling-wave solutions (inden- 
tified when being equals up to translation) at H^, is 



where the existence of such a is provided by the analysis of ([3]) and Z is 
a linear operator from R'^ x M to defined by 



+ P^chiH. I (/3,7)GkerZ} 



Note that on ker Z 



dL,(/3,7) = 7, dc.(/?,7) = A , 

dg.(/5,7) = , di/.(/?,7) = PidcH. + P2hl + /3sd^H, 



41 



and 



Fix (A, u) e and define ^(A, u) a linear operator from 



to 



by 



Ti(/3,7) = -Afch + ^^.'/?i, 

T2(/3, 7) = A [13, < d,H^ >L. +(32 < hi +/?3 < d^H^ >lJ + i^K [M,/?i - 13^] . 

Operator r(A, u) coincide witli (Xdk^, + i/Zc^^dc^,, AdM^, + i'ki,{ki,dCi, — dq^j) on 
kerZ, since as in subsection 14.11 we fiave imposed M^, = H^,(0). 
Now, up to some non-zero F, F', 



D{X,iy) = Fdet(T(A,z/)|kerz) 
and det((r(A, u), Z)) equals 



F'det((r(A,//),Z)) 



ukl 

A < dcH^ >L, +vKM^ \<hl>L, A < ck;H^ >l, 

[dcK] [hi'] [d^Hl] 



-k^u 



-Xkl 




and may be reduced to 
which easily compares with the main part of the Evans function -E'(A, e'^). 



A2 < dcH. >L, +\uk,M, \<hl>L, < dqH, >l, 
\[dMA+^H'M [/iO] [d^H.] 

A [d^l] + vH'!{Q) [/i"'] [d,H':\ 
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